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Mixed-valence phenomena occurring in the "black" (B) and "gold" (G) phases of Smi_i^Ya:S have 
been studied by x-ray diflraction, x-ray absorption spectroscopy, and inelastic neutron scattering. 
Lattice-constant and phonon-dispersion results confirm that the valence instability occurs already 
inside the B phase. On the other hand, pronounced temperature anomalies in the thermal expansion 
a{T), as well as in the Sm mean-square displacements denote the onset of the B-G transition for 
the compositions x = 0.33 and 0.45. It is argued that these anomalies primarily denote an effect 
of electron-phonon coupling. The magnetic spectral response, measured on both powder and single 
crystals, is dominated by the Sm'^"'" spin-orbit component close to 36 meV. A strongly overdamped 
Sm^+ contribution appears only for x > 0.33 near room-temperature. The quasielastic signal is 
strongly suppressed below 70 K, reflecting the formation of the singlet mixed-valence ground state. 
Quite remarkably, the signal around 36 meV is found, from the single- crystal spectra, to arise from 
two distinct, dispersive, interacting branches. The lower peak, confirmed to exist from x = 0.17 
to a; = 0.33 at least, is tentatively ascribed to an excitation specific to the mixed-valence regime, 
reminiscent of the "exciton" peak reported previously for SmBe. 
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I. INTRODUCTION 

Kondo insulators (KIs) form a class of /- (or, for a few 
of them, d-) electron systems which, near room temper- 
ature (RT), typically behave as poor metals, with con- 
duction electrons scattering incoherently off the magnetic 
centers, very much in the same way as in metallic heavy- 
fermions. But contrary to what happens in the latter 
materials, this regime of incoherent Kondo-like scatter- 
ing evolves gradually upon cooling to a semiconducting 
state, due to the opening of a narrow energy gap in the 
electronic density of states below a characteristic temper- 
ature of the order of 50-100 K. This peculiar behavior 
was taken to indicate a strong renormalization of quasi- 
particle states due to electron-electron correlations. For 
this reason, there has been considerable debate^ as to 
whether KIs can be classified as a band semiconductors 
or as Mott insulators, and to what extent their properties 
can be encompassed within the scope of existing theories, 
in particular the periodic Kondo or Anderson models. 
The spin dynamics in these systems is of primary inter- 
est because the latter models imply a direct interplay 
between charge and magnetic degrees of freedom and, in 
some cases, can make quantitative predictions as to, for 
instance, the ratio between a "spin gap" observed in the 
magnetic spectral response and the "charge" gap in the 
electronic band structure. Inelastic neutron scattering 
(INS) studies have now been reported for a number of 



KI materials. The results suggest that one should dis- 
tinguish between systems based on magnetic ions with 
only one unpaired electron (Ce) or hole (Yb) on the 4/ 
shell, and those with a more complex configuration such 
as Sm or Eu. In the following, we will focus on case 
of samarium, for which detailed experiments performed 
on low-absorption isotopically enriched single crystals of 
mixed-valence (MV) SmBg have revealed a number of 
unexpected features;^ In particular, the magnetic re- 
sponse measured at T = 2 K shows, in addition to two 
broad peaks reminiscent of the lower single- ion spin-orbit 
transitions in the Sm^+ (4/^) and Sm^+ (4/^) electron 
configurations, respectively, a very narrow excitation at 
low energy (tiuj « 14 meV) with unusually steep Q de- 
pendence. The latter peak is a genuine feature of the KI 
state, and is strongly suppressed upon heating to only 50 
K. These properties have been explainedii^^ by the for- 
mation of a quantum-mechanically mixed local ground 
state at each Sm site, due to the hybridization of the 4/ 
electrons with the p orbitals of the nearest boron atoms. 
In such an approach, the MV ground state has an es- 
sentially local character and can be regarded as a sort 
of "magnetic exciton" 1^ It can be contrasted with the 
itinerant approach underlying the well-known "coherent 
hybridization gap" picture.'*'^ 

Considering that, in Sm-based materials, both fluc- 
tuating configurations Sm^+ and Sm'^+ have a nonzero 
number of electrons (holes) in the / shell, one may ex- 
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pect exchange interactions between rare-earth (RE) ions 
to play a more important role in the formation of the MV 
ground state than in Ce or Yb. Such effects, however, 
have not been found in SmBg, possibly because of the 
marked local character of the MV state in this system. 
The question of the influence of intersite exchange on the 
ground state of MV semiconductors in general has re- 
mained poorly documented from the experimental view- 
point, with a few exceptions, such as TmSe>i2iii Among 
the Sm compounds, one promising candidate for study- 
ing this problem is SmS. SmS is an archetype unstable- 
valence compound which undergoes a first-order transi- 
tion at room temperature from a divalent semiconduct- 
ing state at P = to a homogeneous MV state above a 
critical pressure of « 0.65 GPa. This transition is ac- 
companied by a volume coUapes, and by a change in 
the sample color from black (B) to gold (G) reflecting 
the metallic conductivity in the MV phase. Penney and 
Holtzbergi?- have shown that the G phase can also be 
obtained at ambient pressure and RT by substituting of 
Y on the Sm site. In the Smi_j.Yj.S crystals they mea- 
sured, the transition was found to take place when the 
yttrium concentration exceeds about 15%, and the lat- 
tice constant falls below ~ 5.70 A. This critical value, 
ac, slightly decreases at low temperature (« 5.67 A at 
T = 80 K). Interestingly, the concentration dependence 
of the lattice constant indicates that, upon substitution 
of trivalent Y, the Sm valence starts to deviate from 2+ 
before the critical concentration, Xc, is reached. This is 
in strong contrast with pure SmS in which Sm remains 
divalent all the way up to the critical pressure. 

In i?-phase SmS at ambient pressure, the magnetic ex- 
citation spectrum at low temperature is associated with 
the spin-orbit (SO) transition from the ^Fo singlet to the 
"^Fi triplet, whose free- ion energy is 36 meV. The de- 
generacy of the triplet is not lifted here because of the 
high. Oh , symmetry of the crystal field at the cubic RE 
site. It was first reported by Shapiro et alw^ thirty years 
ago that the energy of this transition is g-dependent as a 
result of exchange interactions between Sm^"^ magnetic 
moments, mediated by the d-band states. This type of 
magnetic excitations in a system of induced magnetic mo- 
ments is known for RE compounds with a singlet crystal- 
field (CP) ground state, where it is termed "paramagnetic 
(CP) exciton" ji4 The dispersive branch observed in Ref. 
[l3 can thus be characterized as a "paramagnetic SO ex- 
citon" . It was shown to correspond almost exactly to the 
idealized singlet-triplet model. In Ref. ITsL the exchange 
parameters for the first three coordination spheres of Sm 
were obtained from the energy dispersion using a Heisen- 
berg model in the mean-field-random phase approxima- 
tion (MP-RPA). 

Several inelastic neutron scattering (INS) studies have 
been devoted to the phonon spectra of Smi_a;Ya;S, reveal- 
ing a remarkable anomaly in the dispersion of longitudi- 
nal acoustic (LA) phonons in the MV G phases of both 
Smo.75Yo.25S under normal conditionsi^ and SmS under 
pressure«i^ This anomaly was traced back to an electron- 



phonon interaction arising from a low-frequency charge- 
fiuctuation mode specific for the MV state in samar- 
ium. The latter interpretation, based on an excitonic 
description of the charge-fluctuation state first suggested 
by Stevens, ^7 was elaborated in Rcfs. ^ and^ and suc- 
cessfully applied to a quantitative treatment of phonon 
anomalies in Smi_a;Ya;S, G-phase SmS, and SmBg." It 
then served as the starting point for a local description 
of the MV state in Smji as discussed at the beginning of 
this Section. 

There has been, to date, relatively little experimental 
work on magnetic excitation spectra in the Smi_a:Ya;S 
series, and a number of inconsistencies between the re- 
sults have not yet been resolved. In the high-pressure 
experiments of McWhan et al.^, polycrystalline SmS 
was studied on both a triple-axis and a time-of-flight 
spectrometer. The authors report the disappearance of 
the magnetic signal from Sm^+ in the MV G phase. 
However, one can note that the measuring conditions 
were not optimal for the observing an excitonlike (in the 
sense of Ref. W) inelastic peak having a steep Q depen- 
dence of its intensity as found in SmBg. Por the solid 
solution Smo.75Yo.25S, powder experiments at ambient 
pressure22iSi have revealed a broad, Q-dependent, signal 
close to the position of the SO transition in both the 
B and the G phases. In Ref. |^ this broadening was 
interpreted as due to a substructure produced by the hy- 
bridization of magnetic (SO) modes with optic phonon 
modes. Also worth mentioning is the attempt made by 
Weber et al?^ to study the temperature evolution of 
quasielastic magnetic scattering in Smo.75Yo.25S. How- 
ever, their polycrystal, time-of-flight, experiment suffered 
from considerable phonon background, and the measure- 
ments at low temperature were hampered by a low in- 
coming neutron energy, which precluded the analysis of 
the energy-loss side of the spectra. This point will be 
further discussed in section Hvl 

In a previous paper we have reported single-crystal 
INS measurements on the B-phase alloy Smo.83Yo.17S. 
Two dispersive magnetic modes, interacting with each 
other, were observed, lending support to the idea that, 
contrary to the case of SmBg, Sm-Sm exchange inter- 
actions play a significant role in the formation of the 
excitation spectra of the MV state in this system. It 
was argued that the two modes might correspond, re- 
spectively, to excitations from the "parent" Sm^+ state 
and from an excitonlike state reminiscent of SmBg. Pur- 
ther increasing the Y concentration is known to result 
in a transition to the G phase with a higher average Sm 
valence. From the preceding considerations, this should 
lead to i) the suppression of the dispersion, n) enhanced 
excitonic properties of the ground state, implying steeper 
Q- and T- dependences of the corresponding peak inten- 
sities, and in) changes in the parameters of the higher 
energy, single-ionlike excitations associated with the two 
parent states Sm^+ and Sm''+. An important feature of 
the excitonic MV ground state is the correlation found 
experimentally in hexaborides^** between the energy of 
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the new magnetic excitation and the average Sm valence. 
It is not obvious, however, that a similarly strong depen- 
dence of the excitation energy should occur in the sul- 
fides, because / electron states from a given Sm site are 
expected to hybridize with the d states of neighboring Sm 
atoms and, therefore, delocalization may be not so strong 
as in SmBg where hybridization is thought to primarily 
involve the p orbitals from boron. 

In the following, we report a detailed INS study of 
Smi_a;Ya;S (0 < X < 0.45) performed on both pow- 
der and single-crystal samples, covering a wide tempera- 
ture range, and including measurements under an applied 
pressure. We also present x-ray diffraction data on the 
temperature dependence of the lattice constant, the ther- 
mal expansion coefficient, and the Sm mean-square dis- 
placements, as well as the temperature evolution of the 
Sm valence obtained by x-ray absorption (XAS) spec- 
troscopy. A comprehensive analysis of these different re- 
sults allows us to make a connection between dynamical 
characteristics of the material and its valence state, and 
provides clear evidence for the major role played by ex- 
change interactions in the properties of the ground state. 



II. EXPERIMENTS 

The materials used in the present experiments were 
prepared from ^^^Sm isotope with 98.6% enrichment, 
containing less than 0.2% highly absorbing ^''^Sm. Sam- 
ples with Y concentrations x = 0.17±0.03 (a), 0.25±0.02 
(6), 0.33 ± 0.03 (c), and 0.45 (d), as well as Lao.75Yo.25S 
(e) and pure SmS (/) were synthesized, by essentially 
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FIG. 1: (Color online) Room-temperature lattice constants 
of the Smi_i,Ya:S solid solutions studied in this work (error 
bars for lattice constants are less than symbol size); marker 
types for each concentration are used consistently in the rest 
of the paper; closed (open) symbols denote samples which are 
visually in the B (G) phase, respectively, at room tempera- 
ture. Previous data from Ref.lTH (small circles) are shown for 
comparison (see text); dashed line: experimental upper limit 
of the lattice constant for the existence of G phase at T = 300 
K. 



TABLE I: Room temperature lattice parameters ao of mea- 
sured Smi_a,Yi,S and Lao.75Yo.25S samples. 

X 0.17 0.25 0.33 0.45 0.25" 

ao (A) 5.967(2) 5.916(2) 5.835(2) 5.692(2) 5.660(2) 5.769(2) 



identical procedures, either at Tohoku University in 
Sendai (a, c) or at A. F. loffe PTI in St. Petersburg (&, d, 
e, /). Samarium metal was first purified by the Bridgman 
method in a sealed tungsten crucible, then prereacted 
with Y, and S in a silica tube under vacuum using a con- 
ventional electric furnace. Single crystals of pure SmS 
and Smi-a-Ya-S solid solutions were grown by the Bridg- 
man method in electron-beam- welded tungsten crucibles, 
then annealed for several days at 1000°C. Powder for neu- 
tron time-of-flight (TOP) measurements was obtained by 
crushing the smaller single-crystal pieces that were not 
suitable for triple-axis experiments. Smo.55Yo.45S and 
the reference compound Lao.75Yo.25S were synthesized 
in polycrystal form because they were intended for TOP 
experiments only. Powder x-ray diffraction confirmed the 
NaCl crystal structure of all materials, as well as the ab- 
sence of detectable impurity phases. The concentrations 
of samples a, 6, and c were obtained by optical spectral 
analysis, those given above for samples d and e are nomi- 
nal values. Lattice parameters obtained by x-ray diffrac- 
tion at room temperature are listed in Table ^ For the 
compositions x = 0.33 and 0.45, ao is smaller that the 
upper limit of 5.70 A for the existence of the G phase de- 
termined by Penney et alp^ This is consistent with the 
gold color exhibited by crystals c and d. On the other 
hand, the samples with x = 0.17 and 0.25 are found to 
be in the B phase (dark-purple color), with lattice con- 
stants in the 5.8-5.9 A range. As shown in Fig. ^ this 
implies a significant shift in the Y critical concentration 
in the present materials as compared to Ref. Il2l which is 
probably due to differences in preparation procedures. In 
the present data, no discrepancy was found between ma- 
terials from two different origins. Furthermore, the limit 
of ao ^ 5.70 A for the existence of the G phase at RT is 
consistent with all existing results and, in the following, 
ao will thus be assumed to be the relevant parameter for 
comparing results from different studies. 

The temperature dependence of the lattice constants 
in the Smi_2;Ya;S solid solutions (a-rf) and the nonmag- 
netic reference compound Lao.75Yo.25S (e) have been de- 
termined by x-ray diffraction (Cu Ka) in the temperature 
range 10 ^ T ^ 300 K. The values were obtained by re- 
fining the positions of several high-angle reflections, and 
the peak intensities of the same reflections were used to 
derive the relative values of the mean-square displace- 
ments at the rare-earth site. 

The Sm valence for the compositions x — 0.33 and 0.45 
was measured by Lni-edge x-ray absorption spectroscopy 
(XAS) on samples c and d. The experiments were per- 
formed on the station of the DORIS-III storage ring 
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at HASYLAB (DESY, Hamburg, Germany), using a Si 
(111) double-crystal monochromator. The energy resolu- 
tion was estimated to be better than 1.5 eV at 6 keV. The 
measurements were performed by first cooling down the 
sample to 20 K, then increasing the temperature stepwise 
to 300 K. The method for determining the valence was 
described in Ref. l25l 



Neutron powder measurements were performed on 
the TOF spectrometer HET (ISIS, Rutherford- Appleton 
Laboratory, UK) . For each concentration, the total avail- 
able amount of material, of about 4 g, was used to op- 
timize counting statistics. The incoming neutron energy 
was Ei = 300 meV, making the energy range of spin-orbit 
excitations for both Sm^+ and Sm^+ accessible with suf- 
ficient resolution (AE — 13 meV at zero energy transfer, 
further improving in the neutron-energy- loss region). 

Single-crystal inelastic neutron scattering experiments 
were carried out on the 2T thermal-beam triple-axis spec- 
trometer (Laboratoire Leon Brillouin, France). The sin- 
gle crystals measured had volumes of 0.15 (a), 0.35 (6), 
0.17 (c), and 0.42 (/) cubic centimeters. The crystal mo- 
saicity, obtained from the full widths at half maximum 
(FWHM) of the neutron rocking curves on the 200 nu- 
clear peak, was less than 1 degree in all cases. The sam- 
ples were cooled from room temperature (RT) down to 
T = 12 K using a closed-cycle refrigerator. Spectra were 
recorded at a fixed final neutron energy oiEf = 30.5 meV 
(14.7 meV for studying the quasielastic signal), using a 
Cu 111, or PG 002 double-focusing monochromator, a 
PG 002 analyzer, and a graphite filter on the scattered 
beam to suppress higher orders. The resulting resolu- 
tion (FWHM) at zero energy transfer was 2 meV for 
Ef = 30.5 meV, and 0.8 meV for Ef = 14.7 meV. The 
magnetic response was measured in the range of momen- 
tum transfer 1.99 ^ Q ^ 4.6 A"'^. The temperature 
dependence of quasielastic magnetic scattering was stud- 
ied near the Brillouin zone boundary, where this signal is 
well separated from acoustic phonons. A wave vector of 
Q = (1.45,1.45,1.45), with components slightly shifted 
from half-integral values, was chosen in order to avoid 
contamination by second order of Bragg scattering. The 
energy dispersion of selected phonon branches was also 
determined. Data for the longitudinal acoustic branch 
along the [111] direction are presented in Section ITlI Al 

Pressure experiments up to P = 0.48 GPa were per- 
formed using a He-gas chamber cooled down in an ILL- 
type helium- flow cryostat. Measurements on a SmS sin- 
gle crystal (/) and on powdered Smo.67Yo.33S (c), were 
done on the triple-axis spectrometer 2T and the TOF 
spectrometer HET, respectively. In the latter case, a 
near doubling of the incoming neutron flux was achieved 
by reducing the frequency of the Fermi chopper by a fac- 
tor of 2, at the cost of a loss in resolution {AE = 22 meV 
at zero energy transfer). 
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FIG. 2: (Color online) Temperature dependences of lattice 
properties in Smi-j^Y^S as obtained from x-ray diffraction 
measurements: (a) lattice constant for x = 0.33 and 0.45; 
also plotted is the valence of the same compounds obtained 
from XAS measurements; (b) thermal expansion coefficient 
for X = 0.25, 0.33, and 0.45; data for Lao.75Yo.25S are plotted 
as open circles; (c) relative mean-square displacements of RE 
ions at T = 10 K for a; = 0.33 (see text) 



III. RESULTS 

A. Lattice properties 

From the temperature dependence of their lattice con- 
stants, samples a and b {x = 0.17 and 0.25) are found to 
be in the B phase in the whole temperature range from 
RT to r = 10 K. The Sm valence in these samples^^ is 
almost temperature independent, v w 2.19 (a) and 2.22 
(6). On the other hand samples c and d {x = 0.33 and 
0.45) are in the G state at RT, but undergo a phase tran- 
sition at T R:! 200 and w 100 K, respectively, to a state 
with a larger lattice constant [Fig. I^a)] . The tempera- 
tures of the B-G transitions quoted in the following were 
defined from the zero of the second derivative of v(T). 
In Ref. 12, this effect was observed at somewhat lower 
Y concentrations, 0.15 ^ x < 0.3 (see above), and the 
low-temperature (T = 80 K) data points were observed 
to lie on the extrapolation of those obtained for samples 
which remain black at all temperatures, suggesting that 
the low-temperature state is very similar, if not identi- 
cal, to the B phase. Whereas the latter observation was 
not quantitatively confirmed in the case of our samples, 
one can see in Fig. ^a) that the electronic phase tran- 
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FIG. 3: (Color online) LA phonon dispersion along the 
[111] direction for Smo.83Yo.17S and Smo.67Yo.33S; dashed 
hne: data for SmS from Ref. W\ 

sition is indeed accompanied by a considerable decrease 
in the Sm valence at low temperature. Henceforth we 
will use the term "B phase" indiscriminately to denote 
both the low-pressure state and the high-pressure, low- 
temperature state.'*" 

The linear thermal expansion coefficients a = 
{\/a){da/dT) of the above two alloys have been calcu- 
lated from the a{T) data [Fig. IHb)]. Their tempera- 
ture dependence exhibits large additional negative con- 
tributions in comparison with the nonmagnetic com- 
pound Lao.75Yo.25S or the i?-phase alloy Smo.75Yo.25S, 
with pronounced minima located close to the temper- 
atures of the B-G phase transition. From the nuclear 
Bragg peak intensities, we have derived the T depen- 
dence of the Sm mean square displacements, {u'^^*{T) = 
{ul^){T)-{ul^){T = 10 K) in sample c (a; = 0.33). This 
dependence [Fig.l^c)] is nonmonotonic with temperature 
and exhibits a pronounced dip at T ~ 200 K, which co- 
incides with the minimum in the thermal expansion co- 
efficient. The connection between these two effects and 
the electronic phase transition is discussed in Section llVl 

The dispersion curves of LA phonons for Smo.83Yo.17S 
(T = 12 K) and Smo.67Yo.33S (T = 12 K and 200 K) 
along the [111] direction are shown in Fig.|3| Anomalies 
similar to those first reported by Mook and Nicklow^^ 
are observed for all compositions x ^ Q. As expected, 
the phonon softening is particularly strong in the crystal 
with the largest yttrium concentration [x = 0.33) and 
at the highest experimental temperature. However, it 
is remarkable that the effect is still present at tempera- 
tures and Y concentrations where the material is unam- 
biguously in the B phase. Upon decreasing temperature 
and/or Y concentration, i.e. as Sm gets closer to divalent, 
the shape of the dispersion changes and the softening be- 
comes less pronounced. In the Smo.83Yo.17S compound, 
where the Sm valence is almost temperature indepen- 
dent, the LA phonon dispersion does not change with 
temperature within experimental accuracy. This further 
indicates that the phonon anomaly is directly connected 
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FIG. 4: (Color online) Low-temperature (T = 7 K) mag- 
netic scattering function in Smo.83Yo.17S, measured with an 
incident neutron energy of 300 meV at a scattering angle 
{26) = 5° (momentum transfer Q = 1.26 A"^ for E = 2,Q 
meV). The elastic peak, fitted using the spectral lineshape 
obtained from a vanadium standard, has been subtracted out. 
Inset; integrated intensities of the inelastic peaks ai T — 7 
K, normalized to the Sm concentration x = (circle), 0.17 
(diamond), 0.25 (square), 0.33 (triangle), and 0.45 (inverted 
triangle) , plotted as a function of the Sm valence; dashed line 
represents the dependence expected for a weighted superpo- 
sition of magnetic scattering from Sm^"^ and Sm'^'^. 



with the valence instability. 



B. Magnetic excitation spectra in powder samples 

In order to extract the magnetic contribution from the 
experimental TOF data, one first needs to subtract out 
the phonon background, which can be estimated from 
data obtained in a separate measurement of Lao.75Yo.25S 
(sarnple e). The correction procedure was described in 
Ref. l28l and validated by a number of subsequently stud- 
ies. The low-temperature (T = 7 K) magnetic response 
is quite similar for all Y concentrations. Results for 
a; = 0.17 are displayed in Fig. ^ over a wide range of en- 
ergy transfers —20 ^ ii^ ^ 160 meV. The spectrum essen- 
tially consists of a single, relatively broad, peak at about 
34 meV, which can be ascribed to the '^Fo Fi SO 
transition in Sm^+, in agreement with earlier datai^2iSi 
The integrated intensity of the peak, normalized to the 
concentration of Sm ions, decreases when x increases 
(inset in Fig. QIl as a result of the enhanced trivalent 
character. However, the variation shown in the inset ap- 
pears to deviate from a linear slope (dashed line) and to 
level off at the highest Y concentrations. The width of 
the 34-meV peak is comprised between 16 and 18 meV 
for the different samples, which exceeds the experimen- 
tal resolution (FWHM « 11 meV) by a factor of about 
1.5. This broadening cannot be due to simple relaxation 
processes, since the peak shape is not Lorentzian. It can 
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FIG. 5: (Color online) Room-temperature (RT) magnetic 
scattering function in Smi-^jY^S, measured at a scattering 
angle (2^) — 5° (momentum transfer Q = 1.26 and 3.1 
for incident neutron energies Ei = 36 and 130 meV, respec- 
tively). Upper right: expanded plot on the same energy scale; 
data for x = 0.33 at T = 7 K have been added (closed trian- 
gles); solid line: fit of the RT spectrum for Smo.67Yo.33S using 
three peaks (dashed lines) associated with 
{i — 1, 2) (Gaussian), and Sm^"*" ®i?5/2 — *■ "^^7/2 (Lorentzian) 
intermultiplet transitions. 



result from an internal structure of the peak, as observed 
in Ref. I21L and/or an energy dispersion smeared out by 
Q-space averaging in powder data. This interpretation 
was confirmed by our previous single-crystal study of 
Smo.ssYo.irS,^'^ and is further substantiated by the new 
results to be presented in Section UlI CI 

Surprizingly, in spite of the strong mixed- valence char- 
acter of the alloy systems, magnetic intensity could be de- 
tected, at low temperature, neither in the energy range 
of the Sm'^+ ^^^5/2 ^^^7/2 inter-multiplet transition 
(Sg + ~ 130 meV) nor in the quasielastic region. This re- 
sult is important to understand the nature of the mixed- 
valence ground state in this system, and will be discussed 
in more detail hereafter. 

On increasing temperature, strong changes are ob- 
served in the magnetic spectral response. A represen- 
tative set of spectra is collected in Fig. |S1 Extra inelastic 
peaks at energies up to 90 meV, which occur in a similar 
way for all samples studied, can readily be ascribed to 
transitions from the Sm^+ excited multiplets (J = 1,2). 
Furthermore, in the spectra for x = 0.33 and 0.45, a 
very broad signal (FWHM ~ 50 meV), centered around 
E w 130 meV, appears at the highest experimental tem- 
perature (T w 300 K). Its energy roughly corresponds to 
the first SO transition in Sm^+. Under the same con- 
ditions, a similar signal is not observed in samples with 
lower Y concentrations. This can be explained by the 
fact that the samples with x — 0.33 and 0.45 are the 
most trivalent in the series studied, with properties in- 
dicative of the G phase for temperatures higher than 200 
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FIG. 6: (Color online) Magnetic scattering function in 
Smo.67Yo.33S: (a) at P = 3 MPa (T = 10 K); (b) at P = 0.48 
GPa (T = 90 K); solid lines are Gaussian fits to the spectra; 
dash-dotted line through the low-pressure maximum empha- 
sizes the shift in the peak position. Inset in (b): pressure 
dependence of the integrated inelastic magnetic intensity. 



K and 100 K, respectively. It seems therefore natural 
to observe a fingerprint of the Sm'^"'" configuration in this 
regime. However we found no detectable intensity around 
130 meV in the spectrum for Smo.55Yo.45S measured at 
T = 200 K (G phase), even though the valence is higher 
than in Smo.67Yo.33S at room temperature. The reason 
for this discrepancy is not yet understood. Another in- 
triguing issue is the absence of qualitative difference in 
the Sm^+ component between the B and G phases. 

One possibility to stabilize the G phase at low tempera- 
ture within the present composition range, ^ a; ^ 0.45, 
is to apply an external pressure of a few hundreds of 
megapascal. Sample d {x = 0.33), which is close to the 
critical concentration, was selected for this experiment. 
The maximum pressure at Tmin that could be obtained 
in the He pressure cell was 0.4 GPa, which turned out to 
be insufficient to drive the sample into the G phase. As 
can be seen from the comparison of integrated intensities 
for the 34-meV peak measured under different conditions 
(inset in Fig.|HJ), the B-G transition could be reached by 
heating up to T = 90 K, because the melting of the He 
transmitting medium produced an increase of pressure 
to P = 0.48 GPa. The result was a sizable reduction 
in the peak intensity. Since the normal thermal depop- 
ulation of the Sm^+ ground state at 90 K is estimated 
to be no more than 3 per cent, this reduction must re- 
flect a significant decrease in the Sm valence. Another 
interesting observation is that the peak shifts to lower 
energies by more than 3 meV. A similar, though much 
weaker, shift was also observed at ambient pressure with 
increasing Y concentration. The main effects of pressure 
are thus similar to those achieved by Y doping. However, 
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FIG. 7: (Color online) Inelastic neutron spectra in Smi-jjY^iS 
single crystals for Q = (1.3, 1.3, 1.3) at T = 12 K; lines are 
fits by two Gaussian peaks. The ^-Fb — > ^-Fi excitation in pure 
SmS is shown as open circles in (a) (line is a guide to the eye). 



the ^i?5/2 SO transition could not be observed 

under pressure, the reason being probably its weak in- 
tensity and large linewidth, as well as the rather high 
background from the pressure cell. 



C. Single-crystal study 

Three different single-crystals with compositions x — 
0.17 (a), 0.25 (&), and 0.33 (c) were studied in the 
triple-axis experiments. The general evolution of the 
magnetic^^ spectral response in the range of the Sm^+ 
^Fq Fi transition is summarized in Fig. [7| The most 
striking feature, already reported in Ref.|23jfor x — 0.17, 
is the double-peak structure, which is now confirmed to 
exist in the entire MV S-phase regime. One can see 
in Fig. Ul that, for equivalent measuring conditions, in- 
creasing the Y concentration results in a transfer of spec- 
tral weight from the high-energy to the low-energy mode. 
The dispersion curves of the two magnetic branches along 
the main cubic symmetry directions are plotted in Figs. 
|S1 El and ^1 together with the Q dependences of the 
mode intensities. For each concentration, one notes the 
similarity of the energy dispersions for the two branches, 
and the exchange of intensity between them as a function 
of the reduced wave vector q, taking place along all three 
main symmetry directions. Increasing Y concentration 
results in i) a reduction of the energy dispersions for both 
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FIG. 8: (Color online) Wave-vector dependence of the en- 
ergy (upper frames) , and relative intensities (lower frames) , of 
the two magnetic excitations in Smo.83Yo.17S measured along 
main symmetry directions at T = 12 K; lines represent the 
model MF-RPA calculation (see text) 
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FIG. 9: (Color online) Same as Fig.Hlfor Smo.75Yo.25S. 
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FIG. 10: (Color online) Same as Fig.|Hlfor Smo.67Yo.33S. 
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FIG. 11: (Color online) Q dependence of the total integrated 
intensity of the two magnetic excitations in Smi_i:Ya;S alloys 
at T = 17 K along the [111] direction; solid line: calculated 
form-factor for the Fq — > transition; F denotes the Bril- 
louin zone center, L the zone boundary point Q — (|, |, |). 
Intensities have been normalized at the highest experimen- 
tal Q value to the calculated Sm^^ form factor. Inset: par- 
tial integrated intensities of the two magnetic excitations for 
X = 0.33, normalized at Q = 2.66 A~^; open triangles: lower- 
energy peak; closed triangles: higher-energy peak; solid line: 
calculated form factor for the Fo — > Fi transition; wiggling 
lines are guides to the eye denoting the oscillating behavior 
of partial intensities (see text). 
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FIG. 12: (Color online) Quasielastic neutron scattering spec- 
tra for Q = (1.45, 1.45, 1.45) in single-crystal Smo.67Yo.33S at 
T = 16 and 70 K. The spectrum for T = 70 K was fitted 
by a quasielastic Lorentzian lineshape with F ~ 6 meV; solid 
line: fit of the elastic peak at T = 70 K. Inset: tempera- 
ture dependence of the integrated intensity of the quasielastic 
signal; line: calculated temperature dependence from ther- 
mal population effects on the ^-Fq, and F2 multiplets 
for ionic Sm^^, taking into account the T- induced valence 
change; experimental values have been normalized to the cal- 
culated intensities at T = 200 K. 



modes, and a comparable change in their general behav- 
ior; ii) the transfer of part of the total intensity from the 
higher to the lower excitation; iii) an increase in the av- 
erage splitting between the two branches (compare upper 
frames in Figs.l51and ll0|l . The analogy between the dis- 
persions and X dependences of the two modes points to 
a common origin; the fact that they exchange intensity 
further suggests that they may have the same symmetry. 

The integrated intensity associated with the magnetic 
double-peak structure is plotted in Fig. ^] as a function 
of the momentum transfer Q for different yttrium con- 
centrations, X = 0.17, 0.25, and 0.33. The total intensity 
of the two modes decreases with increasing momentum 
transfer, confirming that their origin is magnetic. The 
existence of a superimposed modulation, also observed 
in pure SmS, with a period approximately matching the 
dimension of the Brillouin zone, indicates that they are 
not single-site. The overall Q dependence appears to be 
steeper than the calculated single-ion form factor associ- 
ated with the ^Fo — > Fi excitation in Sm^+, depicted 
by the solid line in the figure. These effects are most pro- 
nounced for X = 0.33. In the inset, the Q dependences of 
the partial intensities for each magnetic mode have been 
plotted separately for the latter composition. They ex- 
hibit a quite remarkable oscillating behavior, again with 
a period equal to the size of Brillouin zone in the direction 
considered. The intensity of the low-energy component, 
averaged over one period, falls off appreciably faster with 
increasing Q than the single-ion form factor, and makes 
the dominant contribution to the observed decrease of 



the total intensity. This gives us a hint that the elec- 
tron states involved in this excitation may have a spatial 
extension larger that the original 4/ shell. 

Another important point concerns the temperature 
evolution of the quasielastic magnetic signal. The most 
detailed data set has been obtained for the composition 
X = 0.33 (Fig. I12|l . In this compound, quasielastic scat- 
tering cannot be detected in the spectrum measured at 
T = 16 K, but a Lorentzian peak, centered aX E — Q, is 
already recovered at T = 70 K. Its spectral weight then 
grows steadily up to room temperature, as shown in the 
inset, while its width increases from ~ 6 meV at 70 K to 
~ 12 meV at 200 K. It is quite notable in the plot that 
the appearance of the quasielastic signal takes place at 
a temperature much lower than expected from the nor- 
mal thermal population of Sm^+ excited SO multiplets 
(shown as a solid line), suggesting that some renormal- 
ization mechanism, associated with the valence fluctua- 
tion, modifies the low-energy electron states within this 
temperature range. 

Finally, we have studied the effect of a moderate hy- 
drostatic pressure on the magnetic excitation spectrum 
in pure SmS. The dispersion curves of the SO excitation 
measured at ambient pressure and for T = 15 K are plot- 
ted in Fig.Hnifor the directions [100] and [111]. The [100] 
data are in excellent agreement with the earlier results 
of Shapiro et al.m^ The dispersion along [111], which was 
not studied in the latter work, can be very well repro- 
duced using the values of exchange parameters obtained 
in Ref. lUfrom a refinement of the [100] and [110] data. 
The energy width of the experimental peak is found to be 
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FIG. 13: (Color online) Energy dispersion of the magnetic SO 
excitation in SmS along two main symmetry directions [100] 
and [111]; open circles: P = 0, T = 15 K; triangles: P = 0.42 
GPa, T = 80 K; closed circles: ambient pressure data at 
T = 77 K from Ref. 1^ dashed lines represent the calculated 
energy dispersions using exchange parameters from Table HTl 

resolution limited. At the maximum pressure of 0.42 GPa 
accessible with the He-gas pressure cell, SmS is still in the 
B phase. As could be expected, no dramatic change in 
the magnetic excitation spectrum occurs with respect to 
P — Q since, in the pure compound, the divalent state 
is retained up to the first-order B-G transition. How- 
ever, measurements performed along the [100] direction 
at T = 80 K (where only the SO ground state is signifi- 
cantly populated) indicate a sizable enhancement of the 
dispersion (Fig.^l left frame). 



IV. DISCUSSION 

A. Lattice properties and valence 

The XAS data reported in Fig.|2fa) for samples c and d 
and in Ref. 26 for the lower Y concentrations clearly indi- 
cate that the Sm valence is nonintegral not only in the G 
phase (samples c at T > 200 K and d at T > 100 K), but 
also in the B phase (samples c and d at low temperature, 
and a and h throughout). This result is in line with the 
conclusion reached by Tao and Holtzberg in their early 
work on Smi-^^Y^^S solid solutions:-" Whereas the nature 
of the valence mixing occurring in the B phase remains 
an open question, the softening of the LA phonon disper- 
sion observed in the present study (Fig.OJ lends support 
to the contention that this mixing is homogeneous in na- 
ture. 

Another interesting issue is the strong lattice expan- 
sion that takes place on cooling in the x = 0.33 and 0.45 
alloys (Fig. I3a)), giving rise to a pronounced negative 
peak in the T dependence of the linear thermal expansion 
coefficient a (Fig. I^b)). These effects indicate a strong 
response of the lattice to the valence transition, with the 
position of the minimum approximately corresponding 
to the transition temperature for each composition. On 



the other hand, despite the significant deviation of their 
valence from 2+, the i3-phase alloys exhibit a normal, 
positive, thermal expansion coefficient, with a T depen- 
dence similar to that in Lao.75Yo.25S. Qualitatively, the 
observed behavior can be explained by the larger ionic 
radius of Sm^+, whose relative fraction increases in the 
average at low temperature. Using the values of the va- 
lence plotted in Fig. |2Ia) , a simple linear interpolation 
("Vegard's law") predicts a volume increase by 1.4 % 
and 1.5 % for x = 0.33 and 0.45, respectively, between 
300 K and 10 K, to be compared with experimental val- 
ues of 2.3 % and 1.3 %. Aside from possible uncertain- 
ties in the above estimates, the difference might point 
to some additional mechanism contributing to the ther- 
mal expansion. In the case of pure SmS, high-pressure 
experiments^i^ have shown the existence, for pressures 
comprised between 0.4 and 1.2 GPa, of a negative peak 
in a{T) whose magnitude is comparable to that found in 
the present work. It can be noted, however, that, unlike 
the lattice constant curves displayed in Fig. |2Ia), those 
reported for SmS under pressure show a normal decrease 
on cooling below room temperature followed by a mini- 
mum around ~ 100 K. As a result, the anomalous depen- 
dence, with negative values of a, develops only at lower 
temperaturesi^ In Refs.l3llandl32. the results were fitted 
to Schottky anomalies, assuming a two-level state with an 
energy gap decreasing gradually from 80-100 K to ~ 30 
K over the experimental pressure range. It is presently 
not clear to what extent such a picture would be appli- 
cable to the Sm2;Yi_2;S alloys because, unlike pure SmS 
at high pressure, these systems exhibit a clear valence 
transition as a function of temperature whose effect on 
the lattice constant certainly cannot be dismissed. 

The present results can be compared with those ob- 
tained previously on SmBg and its alloys with La, which 
also exhibit large negative peaks in their thermal expan- 
sion coefficients. By simultaneously analyzing the tem- 
perature dependences of the heat capacity, Cp, and the 
thermal expansion, one can identify three main terms 
(lattice excitations, /-electron magnetic excitations, gap 
in the electronic density of states), each one being as- 
signed a particular Griineisen coefficient^^ and contribut- 
ing independently (in first approximation) to both Cp 
and a. In pure SmBg, it turns out that the first two terms 
dominate the thermal expansion respectively above and 
below 150 K, whereas the third one is significant only be- 
low 50 K. The important point is that the amplitude of 
the minimum never exceeds — 7xlO^®K^^,as compared 
to —5.5 X 10~^ in the present compounds. This dif- 
ference by almost an order of magnitude reinforces the 
assumption that the conversion of the Sm valence state 
from 3+ to 2+ is the leading mechanism in Smj;Yi_2,S 
alloys undergoing a valence transition as a function of 
temperature, with other thermodynamic effects associ- 
ated with either electronic or lattice degrees of freedom 
making only minor contributions. This fact should be 
kept in mind when comparing sulfide and hexaboride MV 
Sm compounds. 
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As to the Sm mean square displacement {ug^}*(T), 
its non-monotonic variation in Smo.67Yo.33S, with a min- 
imum occurring at about the valence transition tempera- 
ture (Fig-I2jc)), seems to be the consequence of a change 
in the inter-atomic interaction potential due to the B-G 
transition. Such a minimum was not observed, down to 
100 K, in an earlier study by Dernier et al.'^^ on an alloy 
with X — 0.3, and (ugjjj)*(r) was then found to increase 
monotonically with temperature between 100 and 300 K. 
However, in view of the lattice constant data for the same 
alloy reported in Ref. I33L its composition seems to corre- 
spond, in our case, to a higher Y content, close to that of 
sample d. If so, the minimum would have been expected 
to occur near 100 K or below, and thus might have been 
missed in the previous measurements. Let us mention 
for completeness that the sulphur mean square displace- 
ment (Mg}*(r) could not be determined in the present 
experiments, but was found in Ref. "3^ to go through a 
minimum around 230 K. The interpretation proposed by 
Dernier et al. was based on a breathing-shell model in 
which Sm valence fluctuations couple to phonon modes 
having a component of an isotropic deformation of the 
nearest-neighbor S octahedron. It was predicted that 
the only effect giving rise to significant spectral weight 
would consist in a softening of LO modes. The present 
results suggest that the situation is indeed more compli- 
cated and that acoustic modes also may be significantly 
affected. This may provide a basis for further theoretical 
work on the lattice dynamics of MV Sm compounds. 



to T = 16 K, we were able to confirm the complete sup- 
pression of the quasielastic component at low tempera- 
ture. Furthermore, the reappearance of a sizable signal 
already at 70 K cannot be explained by a simple "su- 
perposition" picture, in which Sm^+ and Sm'^+ configu- 
rations would contribute independently. Indeed, there is 
no scattering channel available for quasielastic scattering 
in the J = ground state of Sm^+, and the thermal pop- 
ulation of the excited J = 1 multiplet at T = 70 K is too 
small to explain the intensity observed experimentally. 
In the data of Ref. 0, measured in absolute units, it 
could already be seen that the signal measured exceeded 
the scattering cross section expected from Sm^+ alone. 
The J — 5/2 multiplet of Sm'^"'" could in principle pro- 
vide an additional contribution, but the relatively small 
linewidth observed in Fig. 1121 seems incompatible with 
the damping of the ~* ^-^7/2 transition, which is 

strong enough to completely wipe out the corresponding 
130-meV peak at T = 70 K. Furthermore, if this were 
the case, one would expect the signal existing at 70 K to 
remain visible down to 16 K, which is obviously not true 
experimentally. This leeds us to envision a spin-gaplike 
picture, similar to that proposed previously for SmBg)^ 
in which the ground state is a many-body singlet, differ- 
ent from the single-ion singlet of Sm^+, and forming 
at temperatures below '--^ 50 K. Quasielastic scattering 
then gradually reappears with increasing temperature as 
a result of low-lying excitations from this ground state. 



B. Magnetic dynamics 

1. Quasielastic scattering 

From the neutron scattering results presented in the 
preceding section, the magnetic spectral response of the 
Smi_a;Ya;S alloys was shown to consists of three main 
components: i) a quasielastic signal which exists only 
for temperatures higher that ~ 50 K; m) a doublet 
of peaks in the energy range expected for the Sm^+ 
Fq — > Fi inter-multiplet transition (with extra peaks 
attributable to transitions from excited multiplets ap- 
pearing at increasing temperature); Hi) only in the G 
phase and at high enough temperature: a very broad 
signal around 130 meV, which can be traced back to the 
Sm3+ ^i75/2 transition. 

The temperature dependence of the quasielastic signal 
for X — 0.33, summarized in Fig. 1121 corresponds to a 
temperature range below the transition, which occurs at 
~ 160-180 K. The system is thus expected to be predom- 
inantly in the B phase. The variation of the intensity 
shown in the inset turns out to be in good agreement 
with the TOF resuhs (50 K < T < 250 K) previously 
reported by Weber et al,-^ for a sample with x = 0.25. 
In the latter study, the B-G transition presumably took 
place around 150 K but no clear effect was noted on the 
neutron data. By extending the range of measurements 



2. Inter-multiplet transitions 

One remarkable feature of the magnetic response of 
the Smi_2:Y2:S solid solutions is that, whereas we have 
found the low-energy part of the magnetic response to 
be deeply altered as a result of valence mixing, the inter- 
multiplet transitions corresponding to Sm^+ are observed 
at about the same energies (w 35 meV for '^Fq ^Fi) 
as in divalent compounds, notably SmS itself, and with 
the same type of energy dispersion as in the latter com- 
pound. There is, however, one major difference, which 
is the appearance of a second excitation located about 5 
meV lower in energy than the original peak. Already in 
the early TOF experiments of HoUand-Moritz et al^ on 
Smo.75Yo.25S powder, indication was found for a struc- 
ture in the inelastic neutron peak, implying the existence 
of two or more distinct modes. This effect was tenta- 
tively ascribed to a magnetoelastic effect, assuming a 
coupling to occur between the dispersive Sm^"*" "spin- 
orbit exciton" branch and an optic phonon branch of 
the same symmetry. Evidence against this interpretation 
now comes from the observation that this doublet peak 
structure does occur in the compound with x — 0.17 even 
though, for this composition, the frequencies of the op- 
tic phonons lie well below the spin-orbit exciton band, 
making the suggested electron-phonon process basically 
ineffective. Furthermore, the Q-dependence of the inten- 
sities for the two peaks (taken either separately or as a 
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whole) globally follows that expected for a pure magnetic 
process (single-ion form factor), ruling out a significant 
admixture of phonon character. 

Another possible mechanism which might account for 
the double peak structure is a crystal-field (CF) split- 
ting of the J = 1 Sm'^+ multiplet«2i Formally, the local 
symmetry lowering expected to result from the substitu- 
tion of divalent Sm ions by Y'^^ on the first rare-earth 
coordination sphere can produce non-cubic components 
in the CF potential, which ought to result in a splitting 
of the ^i^i excited triplet state. For instance, a tetrag- 
onal configuration of the charge environment will lead 
to a splitting of this triplet into two substates: Jz = 
and Jz = ±1. However, this mechanism cannot explain 
the two-peak profile observed in the magnetic excitation 
spectra. The point is that, for all Y concentrations stud- 
ied, there is no dominant configuration of the first rare- 
earth coordination sphere around each Sm^"*". Even for 
X = 0.17, the probabilities to find configurations with 
one, two, and three Y ions among the 12 nearest neigh- 
bors are 26, 30, and 20 per cent, respectively. Each of 
these configurations, which in total represent about three 
quarters of the sites and should thus dominate the mag- 
netic response, gives rise to different local distortions of 
the CF potential depending on how the Y ions are dis- 
tributed amongst the 12 available sites. From such a 
picture, the occurrence of two well defined peaks seems 
very unlikely. This argument is further supported by 
previous neutron measurements— on the inhomogeneous 
MV compound Sm3Te4, in which inter-atomic distances 
are close to those in SmS. In that compound, the ratio 
of Sm^+ to Sm^+ is exactly one half, and the probabil- 
ities of finding 2-\- and 3+ ions on the first rare-earth 
coordination sphere of Sm are comparable. Nevertheless, 
the peaks in the neutron spectra are well defined and 
corresponds to the CF scheme corresponding to the orig- 
inal, undistorted, local symmetry of the rare-earth site 
with chalcogen atoms as the first nearest neighbors. The 
only indication of an effect of the CF distortion due to 
charge inhomogeneity was possibly a weak broadening of 
the peaks in the neutron spectra. Therefore, it is rea- 
sonable to assume that, in Smi_j;Yj;S too, the effect of 
non-cubic components of the CF potential is relatively 
unimportant, and can at most produce some broadening 
in existing SO excitations. 

Other peculiarities of the extra mode, in particular the 
steep Q dependence of its intensity displayed in Fig.lTTlor 
the energy dispersion, indicative of a coherent effect and 
thus a priori difficult to reconcile with a random local en- 
vironment, are also against a CF origin. But at the same 
time, they might provide some clues as to what could 
be a more plausible mechanism to account for the neu- 
tron data. In Rcf. 23, a simple phenomenological model 
was proposed, based on the assumption of two magnetic 
modes hybridizing with each other, whose dispersions can 
be represented in the formalism of a singlet-triplet mean- 
field-random-phase-approximation (MF-RPA) model. In 
the context of the "local-bound-state" picture of Ref. i^, 



these modes are assigned respectively to the "parent" 
Sm^+ ^Fq Fi spin-orbit transition, and to a renor- 
malized excitation (J* = ^ J* = 1) associated with 
the quantum-mechanical MV state. Other physical ori- 
gins of the second mode are in principle compatible with 
the model, provided it retains the right symmetry to hy- 
bridize with the spin-orbit exciton, but the present as- 
sumption provides a natural explanation for the steeper 
Q dependence of the intensity because the local bound 
state is less localized than ionic 4/ states. Let us note, 
however, that in contrast to the case of SmBg, the in- 
tensities of both branches do not follow the variation of 
the (local or extended) single-site form factor because of 
exchange couplin g an d mode-mode interaction. 

Following Ref. 13, the dispersions of the two magnetic 
modes {X — a, (3 for the "parent" and renormalized spin- 
orbit excitations, respectively) in the noninteracting limit 
are given by 



1 - 2MlR{T) 



^(q) 
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where is the squared matrix element of the tran- 
sition, J(q) the Fourier transform of the Sm-Sm in- 
direct exchange interaction, and R{T) a temperature 
factor. The two modes interact with each other ow- 
ing to f — d mixing or partial delocalization of / elec- 
tron, which is also responsible for the existence of the 
lower mode. Therefore mode-mode interaction can be 
regarded as the consequence of two effects: MV state 
formation (on-site effect) and / — / exchange interac- 
tion (inter-site effect). The mixing interaction between 
the two modes is thus written as the sum of a (domi- 
nant) on-site term Vq and a (smaller) hybridization term 
^i(q) — J2^NN exp[ici{rj — ro)] due to the hybridiza- 
j 

tion (Vnn) of a mode from a given Sm at site Tq with 
the other mode from a second Sm at site Vj belonging to 
its first coordination sphere. The above interaction can 
be treated in standard second-order perturbation theory, 
yielding the energies and intensities for both hybridized 
modes. 

In Ref. 23 the experimental dispersion curves for 
Smo.83Yo.17S for the the two modes were fitted using 
only two adjustable parameters Vq and Vi. Interestingly, 
it was found that using the parameters derived from the 
fit of the energy dispersions, the above phenomenologi- 
cal model also qualitatively explains the transfer of in- 
tensity between the two modes as observed experimen- 
tally. In fact, the exchange interaction alone is sufh- 
cient to roughly reproduce the shape of the dispersion 
curves but the hybridization term is necessary to ac- 
count for the variation of the intensities. In the calcula- 
tion of J(q), the exchange constants up to third nearest 
neighbors were considered following Ref. Il3l and their 
values were assumed to be the same for both modes, 
equal to those previously determined in undoped SmS. 
Aq was taken equal to 36.2 meV, the single-ion energy 
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TABLE II: Parameters of the model used to describe the dispersive magnetic modes in Smi-^Y^iS. The indices a and (3 refer 
to the two different excitations observed experimentally (see text). 



X 


Sm 




Number of 




Energy 


Exchange constants 


Mixing parameters 


Relative intensities 




valence 


1st neighbors 


(meV) 




(meV) 






(meV) 


of bare excitations 






ni 


n2 


ns 


Ed 




Ji 




Js 




Vi 






0" 


2.00 


12 


6 


24 


36.2 




0.043 


0.025 


-0.003 













2.00 


12 


6 


24 


36.2 




0.043 


0.020 


-0.003 










0.17 


2.20 


10 


4.95 


20 


36.2 


32.5 


0.050 


0.020 


-0.003 


-0.60 


0.25 


0.58 


0.42 


0.25 


2.22 


9 


4.50 


18 


36.2 


31.5 


0.054 


0.020 


-0.006 


-0.60 


0.25 


0.55 


0.45 


0.33 


2.34 


8 


4 


16 


36.2 


30.5 


0.056 


0.025 


-0.020 


-0.60 


0.25 


0.44 


0.56 



"Ref.ES 



of the Sm2+ ^Fq -> '^Fi transition, and to 32.5 meV, 
from the empirical relation found in the S mi^j^La^^Bg se- 
ries of alloys between the local bound state energy and 
the Sm valence;^^ In order to consistently analyze the 
data from all three samples, the fitting procedure applied 
previouslji^ was modified so as to to take into account 
the general trend of q dependencies for both energies and 
intensities in the entire series. This leads to the fitting 
curves drawn in Figs.lSl- irUI with parameter values listed 
in Table in 

If we first consider the variation of exchange interac- 
tions, it appears that, as the concentration of nonmag- 
netic substituent increases, with an attending decrease 
in the lattice spacing, the antiferromagnetic contribution 
rises substantially (almost by a factor of 7), producing 
a shift of the minimum in the energy dispersion curves 
from g = to a finite q value, as can be seen in the fig- 
ures. It is worth noting that the application of an exter- 
nal pressure to pure semiconducting SmS results in just 
the opposite effect on the exchange parameters, namely 
an enhancement of the ferromagnetic character in the 
dispersion curves. In Fig El the energy at the F point 
(Brillouin zone center) appears to be lowered by approx- 
imately 1 meV at a pressure of 0.42 GPa. This corre- 
sponds to an increase of the Sm-Sm exchange coupling 
by more than 10 per cent. The effect of the mere re- 
duction in the lattice spacing, calculated using published 
SmS compressibility dataj^ is found to be only about 2 
per cent. Therefore the enhancement of the interaction 
between Sm ions (accounted for by the exchange inte- 
grals) could be supposed to take place under the applied 
pressure, in analogy with a Eu chalcogenidesi^ This fact 
indicate some differences existing between effects of phys- 
ical and chemical pressure on the one hand, and between 
inter-ion interaction in normal and mixed valence state 
on the other one. 

Regarding hybridization parameters (normalized in- 
tensities of the bare excitations A\ = Ml/(^M^ + -^^j) 

and mode- mode interaction V{q)), the variation with Y 
doping cannot be uniquely ascribed to one particular pa- 
rameter without making an extra assumption. In the 
present analysis, it was supposed that the main effect 
is a change in the relative weight {Ap/Aa) of the two 
modes, while other parameters (Vq, Vi) were fixed at the 



same values as in Ref. |23 for x = 0.17. The results re- 
ported in Table ITll indicate a significant enhancement of 
the lower "renormalized" excitation. It can also be noted 
that, in Smi-j^Y^jS, the renormalization of the energy 
of the lower peak as valence increases is not so strong 
as for the "exciton" peak in Smi.^^La^Bg. Whereas for 
Smo.83Yo.17S an optimum fit is achieved using the value 
of calculated from the empirical relationship derived 
in hexaborides'^ between the peak energy and the valence, 
the change in this value required to fit the data for larger 
concentrations x = 0.25 and 0.33, is not as large as would 
be expected from the latter relationship. 



The above description of the single-crystal data is in a 
good agreement with the time-of flight data. Specifically, 
it gives us clues to some peculiarities of the magnetic 
spectral response measured on powder, like for instance 
the shift of the centroid of the 34 meV peak produced by 
Y doping or, even more clearly, by the application of pres- 
sure (see FigEJ. This effect can now be traced back to 
the enhancement of the low-energy subpeak, correspond- 
ing to the renormalized spin-orbit transition. The doping 
dependence of the integrated peak intensity (inset in Fig 
0J can also be explained by the MV character of the extra 
magnetic mode. If the existence of this mode were simply 
due to a CF splitting of the original spin-orbit transition, 
this would result in a simple redistribution of the original 
''i^o "^Fi cross section between the two satellites, and 
the total integrated intensity should decrease proportion- 
ally to the deviation from the divalent state. As already 
noted, this is not verified experimentally. This discrep- 
ancy may be solved if one considers that the lower peak, 
which becomes stronger as the valence increases, has a 
different origin, with additional magnetic degrees of free- 
dom getting involved as a result of the mixing with more 
extended electron states. Increase of the temperature is 
expected to suppress excitonic transitions in the spec- 
trum quite fast. Indeed, Y doping at room temperature 
(or temperature increase for the particular Y concentra- 
tion, x ( 0.33) does not result in some visible shift of this 
peak. 
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V. CONCLUSION 

Smi^ajYj^S solid solutions with compositions a; = 0, 
0.17, 0.25, 0.33, and 0.45 have been studied by x-ray 
diffraction and inelastic neutron scattering, both on pow- 
der and single crystals. By relating lattice properties 
(thermal expansion coefficient a{T), mean square dis- 
placement (wgni)) to the Sm valence obtained from XAS 
spectra, it was shown that all samples with x > 0.17 
are in a homogeneous mixed- valence state. All measured 
samples are in the B phase at low temperatures, and 
the G phase is found only at T > 200K {x — 0.33) and 
T > lOOK {x = 0.45). The B-G electronic phase tran- 
sition is accompanied by a change in the inter-atomic 
interaction potential, thereby affecting the temperature 
dependence of a and {ug^). Further evidence for the 
anomalous character of the B phase, even at the lowest 
Y content x = 0.17, is provided by the clear anomaly 
in the LA phonon dispersion curve (previously reported 
only for the G phase). 

For all concentrations studied, the magnetic excitation 
spectrum of powder samples displays a broad peak at 
energy £^ « 34 meV corresponding approximately to the 
Sm^+ ^Fg Fi spin-orbit transition. As Sm becomes 
more trivalent, this peak weakens and, at high enough 
temperatures and for valences exceeding ~ 2.45, a very 
broad signal, reminiscent of the ^i?5/2 — > ^^^7/2 transi- 
tion in Sm'^^, is observed around 130 meV . This strongly 
damped response likely indicates a very high rate of spin 
fluctuations. For T > 70 K, an additional quasielastic 
signal appears. As in SmBg, its temperature dependence 
cannot be explained by simple thermal population ef- 
fects, and can be traced back to the emergence of the 
quantum mechanical mixed-valence ground state below 
approximately this temperature. 

The most remarkable feature of the magnetic spectral 
response of Smi_3;Y3;S in the mixed valence regime was 
established by the single crystal experiments, namely an 
extra excitation occurring on the low-energy side of the 



Sm^+ ^Fq — > Fi spin-orbit transition. The evolution 
with increasing Sm valence of the intensities and energy 
dispersions of both spectral components can be success- 
fully reproduced by a simple phenomenological model 
based on two hybridized modes. The formation of a 
loosely bound /-electron state ("MV exciton"), devel- 
oped in detail for the case of mixed- valence SmBg,^ is 
proposed as the possible physical origin of the novel ex- 
citation. The main difference with SmBg is the q depen- 
dence of the energies and intensities for both mode. This 
effect is ascribed to Sm-Sm exchange interactions via the 
d band, which is simultaneously involved in the process 
of /-electron delocalization due to valence instability. 

In conclusion, the present results shed some light 
on lattice and magnetic aspects of the long-debated 
problem of valence fluctuations in Sm systems. In 
particular, the INS results on magnetic dynamics reveals 
clear similarities, as well as conspicuous differences, 
between the sulfide alloys and the previously studied 
hexaboridc compounds. In this connection, the discovery 
of properties reminiscent of the latter materials in a 
completely new family, Sm fuUeridej.™ opens promising 
possibilities to check our present understanding of these 
phenomena. 
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